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Abstract 

The interface stability against small perturbations of the planar solid-liquid in- 
terface is considered analytically in linear approximation. Following the analytical 
procedure of Trivedi and Kurz (Trivedi R, Kurz W. Acta Metall 1986;34:1663), 
which is advancing the original treatment of morphological stability by Mullins 
and Sekerka (Mullins WW, Sekerka RF. J Appl Phys 1964;35:444) to the case of 
rapid solidification, we extend the model by introducing the local nonequilibrium 
in the solute diffusion field around the interface. A solution to the heat- and mass- 
transport problem around the perturbed interface is given in the presence of the local 
nonequilibrium solute diffusion. Using the developing local nonequilibrium model 
of solidification, the self-consistent analysis of linear morphological stability is pre- 
sented with the attribution to the marginal (neutral) and absolute morphological 
stability of a rapidly moving interface. Special consideration of the interface stabil- 
ity for the cases of solidification in negative and positive thermal gradients is given. 
A quantitative comparison of the model predictions for the absolute morphologi- 
cal stability is presented with regard to experimental results of Hoglund and Aziz 
(Hoglund DE, Aziz MJ. Mat Res Soc Symp Proc 1992;205:325) on critical solute 
concentration for the interface breakdown during rapid solidification of Si-Sn alloys. 
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1 Introduction 



In the solidification of liquids, an initial solid-liquid interface is subject to 
unstable growth which leads to various crystal patterns of cellular, dendritic, 
banded, fractal, etc. morphologies [1,2,3,4,5]. A scheme of changing the crystal 
morphology with the interface velocity can be considered in the example of 
one-phase solidification, i.e. when the liquid transfers into solid without pre- 
cipitation of additional phases. Figure 1 shows schematically the steady-state 
growth morphologies which form in a liquid as a result of morphological insta- 
bility at given interface velocity, V, in single-phase solidification. With a small 
velocity, an initially smooth interface remains planar up to a velocity equal to 
the critical velocity, Vc, defined by the constitutional undercooling [6]. Behind 
Vc, the smooth interface becomes unstable and the interface exhibits a steady 
cellular morphology. By further increasing the velocity, a surface of cells may 
become unstable with the developing of dendritic patterns. At high interface 
velocity, dendritic patterns degenerate with the appearing of rapidly moving 
cells. A demarcation line at V — Va divides the regions between the interface 
instability, V < Va, and the absolute stability, V > Va, where the planar in- 
terface is morphologically stable against small perturbations of its form. This 
demarcation line is usually known as the critical velocity, Va, for absolute sta- 
bility of the planar interface. The sequence of growth morphologies, Fig. 1, is 
well-known from experiments on directional solidification and solidification in 
the undercooled state [7]. It has been demonstrated in computational modeling 
[8] of crystal growth as well. 

The theory of morphological stability was developed first by Mullins and Sek- 
erka, who considered the stability of a spherical particle grown into a su- 
persaturated solution [9], and stability of the planar front during directional 
solidification of a binary liquid [10]. In investigating the evolution of small 
harmonic perturbations of the planar interface [10], they provided a rigorous 
basis of linear morphological (in) stability at low solidification velocity. Partic- 
ularly, Mullins and Sekerka introduced a concept of marginal stability for the 
wavelength of perturbation, which gives neutral stability of the plane when 
the amplitude of perturbation does not change in time. Behind the critical ve- 
locity of the marginal stability, the increasing of the amplitude of perturbation 
in time may lead to cellular or dendritic interfaces. The theory [10] gave rise 
to a great number of investigations of morphological transformations due to 
linear instability of interfaces and nonlinear behaviour of unstable interfaces. 
These are presented in overview by Coriell and McFadden [11] and monograph 
by Davis [12]. 

In its classic form [10], the theory assumes a local equilibrium at the interface, 
which is an excellent approximation for many metallic systems solidifying at 
small interface velocities. At large driving force for the interface advancing, 
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and with the increasing of its velocity, the analysis of Mullins and Sekerka 
can be modified. Trivedi and Kurz [13] extended the analysis of Mullins and 
Sekerka [10] for the case of rapid solidification, and introduced the stability 
functions dependent on the interface velocity. By taking into account the veloc- 
ity dependent coefficient of solute redistribution (partitioning function), they 
developed an analytical model [14] of microstructure formation under direc- 
tional solidification, over the range from low interface velocity up to velocity, 
Va, of absolute morphological stability. 

In addition to the previous analysis of pattern formation and morpholog- 
ical stability of the interface in which the treatment is extended to rapid 
solidification [13] and nonequilibrium effects at the interface [14], the local 
nonequilibrium in bulk phases may play an essential role in rapid solidifica- 
tion. Particularly, a deviation from local equilibrium in solute diffusion may 
act on the rapid advancing of the solid-liquid interface because the interface 
velocity, V, can be of the order or even greater than the solute diffusion speed, 
Vd, in bulk phases. For instance, the diffusion speed can be of the order of [15]: 
Vd ~ 0.1 — 10 m/s. In modern experiments on solidification of undercooled 
droplets the interface velocity approaches [16]: V = 10 — 100 m/s. Therefore, 
the undercooling of liquids is sufficient for detecting solidification with the 
interface velocity comparable to the diffusion speed. 

Considering the process of solute diffusion, Fick's first law is obtained on the 
basis of classic irreversible thermodynamics of Onsager and Prigogine which 
assumes propagation of concentration disturbances with an infinite speed and 
with local equilibrium in the bulk phases [17]. Local equilibrium is character- 
ized by the statistical distribution function, given by the first order term of 
its expansion [18]. However, for a high- velocity solidification front, the time 
for crystallizing of a local volume is comparable to the time for relaxation 
of the diffusion flux to its steady-state value [15]. In this case, local equilib- 
rium is absent in the bulk phases and the solute flux cannot be described by 
the classical Fick's first law. Including the evolution equation for the diffu- 
sion flux, the analysis of Galenko and Sobolev [19] shows that the deviations 
from local equilibrium in phases and at the interface drastically affect both 
the solute diffusion and the interface kinetics. Therefore, in this article we 
consider the linear stability analysis for a rapidly moving interface under local 
nonequilibrium solute diffusion. 

The linear morphological analysis of the interface stability has been developed 
and used to obtain a stable tip of dendrite growing under local nonequilib- 
rium solute diffusion in rapid solidification [20]. A marginal stability criterion 
was used and concluded that at V — Vd the complete transition to diffu- 
sionless solidification may proceed sharply, with the appearing of the break 
point in the kinetic curves "dendrite velocity - undercooling" and "dendrite 
tip radius - undercooling". In addition to this, Lee et al. [21] performed the 



3 



linear stability analysis in rapid directional solidification using the model of 
Galenko and Sobolev [19]. As they showed [21], the effect of the local nonequi- 
librium in solute diffusion postpones the onset of the cellular instability in 
better agreement with experimental data, in comparison with the model pre- 
dictions in which the local nonequilibrium only at the interface is considered. 
However, the analysis of Lee et al. [21] was not self-consistent (in the analysis, 
the authors used an expression for the slope of kinetic liquidus obtained from 
the local equilibrium thermodynamics). As it has been noted in Ref. [20], the 
predictions of rapid solidification of alloys can be satisfactory compared with 
experimental data only on the basis of self-consistent model (i.e., when all 
model functions are taking into account the deviation from local equilibrium 
in the solute diffusion field). Consequently, the first purpose of the present 
article is to analyze the morphological stability of the planar interface on the 
basis of the self-consistent model of local nonequilibrium solidification. 

In the present analysis, special consideration is paid to the stability of the 
interface around the demarcation line, V = Va, below which the interface 
instability occurs and behind of which the absolute morphological stability 
proceeds, Fig. 1. This transition is known from experiment as a transition 
from cellular patterns to the segregation-free patterns [22,23]. The crystal 
microstructure after the transition is shown in Fig. 2. As it has been stated [24], 
the transition from the macroscopically smooth solid-liquid interface to the 
cellular-dendritic microstructures occurs with the decreasing of the interface 
velocity, V, below critical velocity, Va, for absolute stability of the planar 
front. Around the velocity V — Va, the kinetics of crystal growth begin to 
disagree with the predictions of the model in which the local nonequilibrium 
only at the interface is considered (see the results of analysis in Ref. [20]). 
Consequently, the second purpose of this article is a quantitative evaluation 
of the discrepancy between the present local nonequilibrium model and the 
model in which the local nonequilibrium only at the interface is taken. It is 
also in comparison with experimental data on morphological stability of the 
interface in solidifying alloy. 

The article is organized as follows. In Sec. 2 we give a formal description of 
departure from local equilibrium due to solute diffusion and give a set of gov- 
erning equations to analyze the morphological stability of the interface. In Sec. 
3 we analyze the influence of perturbations in fields and at the interface on 
the linear stability of the planar interface. In Sec. 4, an obtained criterion of 
marginal (neutral) stability allows us to analyze the morphological stability of 
the interface for the case of solidification an undercooled melt (with negative 
temperature gradient) and for the directional solidification (with the positive 
temperature gradient). The absolute stability of the planar interface is ana- 
lyzed in Sec. 5. A discussion about expressions for solute trapping and kinetic 
liquidus which define a final form of the function for the absolute stability 
of the planar interface is given in Sec. 6. Also, in this section, we compare 
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the derived function for the absolute interface stability with the available ex- 
perimental results obtained in Si-Sn alloy solidification. Finally, in Sec. 7 we 
present a summary of our conclusions. 



2 Statement of the problem 

We shall consider a dilute binary alloy that undergoes nonisothermal solidifi- 
cation in infinite space. Let us take into account the heat diffusion in phases, 
solute diffusion in the liquid, and one can neglect the solute diffusion in solid. 
The main physical assumption of the present problem is an absence of local 
equilibrium both at the solid-liquid interface and in the solute diffusion field 
around the interface. In this case, the degree of local nonequilibrium is esti- 
mated by the relation of the interface velocity, V, and the diffusion speed, Vd, 
which is a parameter of the process of diffusion and can have different values 
at the interface and in bulk phases. The speed Vd is a maximum speed of prop- 
agation of the diffusion profile in the system and defined as Vd = (D/td) 1 ^ 2 , 
where D is the diffusion coefficient, and td is the time of relaxation of diffu- 
sion flux to its steady-state value. Therefore, we develop the rapid solidification 
model which is taking into account the finiteness of the diffusion speed in the 
system. 

2. 1 Departure from local equilibrium 

If local thermodynamic equilibrium in the bulk is not reached, the connection 
between the vectors of diffusion fluxes, and J, and the driving forces, VTj 
and VC, for the heat and solute diffusion, respectively, have the following 
integral form: 

- relaxation of the heat flux 

q*M) = -f Di(t-t*)VTi(t*,T)dt*, (1) 

- relaxation of the solute diffusion flux 

J(r, t) = - f Dj{t - t*)VC(t*, r)dt*, (2) 

J — oo 

where index i = L or i = S is related to the liquid or solid phases, respectively. 
Tj are the temperatures in the phases, C is the solute concentration in the 
liquid, t is the time, r is the radius-vector of a point in the system, and 
Dnit — t*) are the relaxation functions of the fluxes (R — q or R — j). 
Equations (1) and (2) imply the fact that when the interface moves with a 
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high velocity, local equilibrium in the fields does not occur and the diffusion 
fluxes at a point in the system no longer depend on the instantaneous gradients 
of the temperature and chemical composition, but are also determined by the 
local prehistory of the solidification process. 

Equations (1) and (2) represent general expressions for evolution prehistory of 
the diffusion processes. For the case of heat diffusion, when the heat propagates 
with much higher speed in comparison with the interface velocity, the influence 
of local nonequilibrium in the temperature field on the kinetics of the interface 
advancing is negligible (see the analysis of heat transfer in rapid solidification 
in Ref. [25]). Therefore we specially define the relaxation functions Dr in 
Eqs. (1) and (2) for the important class of dissipative hyperbolic models in 
which they take the following forms 

D\(t-lT)=D\{0)6{t-t*), (3) 

D,.(*-O=^-(0)exp(-^), (4) 

where D l (0) = Ki are the thermal conductivity in the liquid (i = L) and 
solid (i — S), S is the Dirac delta-function, and Dj(0) = D/t d is a value of 
relaxation function for solute diffusion at a moment t = t*. 

Eq. (3) describes an instant relaxation which occurs at a moment t = t* . 
Therefore, one can expect a description of local equilibrium heat transport by 
the function (3) in combination with Eq. (1). In contrast to this, the relax- 
ation function (4) with the flux prehistory (2) leads to local nonequilibrium 
solute diffusion. In usual circumstances, the relaxation time, td, for diffusion 
flux is very small and the relaxation effects are negligible. It is usual to speak 
in this case of "viscous diffusion" which can be accurately described by the 
Fickian local equilibrium approximation. However, equations (2) and (4) de- 
scribe relevant difference of solute diffusion with respect to the classical Fick's 
law. To date, departures from this law are known as the appearance of inter- 
nal effects, couplings of diffusion and viscosity, and longitudinal diffusion (see 
Refs. [26,27,28,29] and references therein). In addition to these appearances, 
Eqs. (2) and (4) can also be applied to the process of the phase transition 
in a strongly nonequilibrium medium. Rapid solidification of a binary melt is 
a good example of such a nonequilibrium phase transition in which the in- 
terface can move with the high velocity comparable with the diffusive speed, 

V ~ V D ~ 0.1 — 10 (m/s) [15]. In this case, the time for crystallizing of local 
bulk is of the order of the relaxation time of the solute diffusion flux [15,19,20], 
characterizing its decay forward of its local equilibrium value. Therefore when 

V ~ Vd, the relaxation interacts with the diffusion process directly, and it is 
necessary to take into account the local prehistory of the solute diffusion, e.g. 
in a form described by Eqs. (2) and (4). 
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Eq. (4) simulates a physically reasonable situation in which exponential decay 
of the diffusion flux occurs in the local bulk of the liquid phase. This equation 
provides the lowest order of approximation of the diffusion flux relaxation. 
Indeed, substituting Eqs. (3) and (4) for Eqs. (I) and (2), respectively, one 
obtains 

q, + KiVTi = 0, (5) 
81 

r D — + 3 + DVC = 0. (6) 

Equation (5) is a well-known Fourier law which is true for infinite thermal 
speed in the system, i.e. the heat diffusion flux is instantly relaxed to its local 
equilibrium value and the effects of local nonequilibrium in the thermal field 
are negligible. Equation (6) can be treated as the simplest generalization of 
the classical Fick's first law J + DVC = that is recovered when td = or 
in stationary situations in which dJ/dt = 0. The evolution equation (6) takes 
into account the relaxation to local equilibrium of the mass flux and is known 
as the Maxwell-Cattaneo equation in the context of heat transport [27,28,29]. 
By taking the relaxation (2) with the exponential law (4), it follows from 
Eq. (6) that the flux J at a point in the system is defined by the evolution of 
the concentration gradient VC(f ,r) during the period t — r D < t* < t, but 
not by the gradient VC(t, r) at the moment t, as in local equilibrium approxi- 
mation. Thus, taking into account the exponential decay for the relaxation of 
the diffusion flux (see Eq. (4)), the simplest evolution equation (6) is obtained 
for the interaction of the relaxation process and solute diffusion. Instead of 
Eq. (6), by taking suitable relations for the functions of D(C), D(dC/dt), or 
D(VC), one can describe more complicated situations for non-Fickian diffu- 
sion in nonequilibrium media. These are described by the evolution equation 
with the higher order time derivatives or couplings of relaxation to non-local 
effects for transient processes 



The heat and solute diffusion are governed by the balance laws 

dT- 

X ^ + V-q, = 0, (7) 

dC 

^ + V.J = 0, (8) 

where Xi are the heat capacities in phases. Substitution of Eqs. (5) and (6) 
into Eqs. (7) and (8), respectively, gives the following system of equations 

^ = a,V%, (9) 
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where Oj are the thermal diffusivities in the liquid (i = L) and solid (i = S). 
Equation (9) is the common partial differential equation of parabolic type for 
the heat transfer which adopts the infinite thermal speed. Equation (10) shows 
that Eqs. (6) and (8) give rise to the partial differential equation of a hyperbolic 
type for the solute concentration, which is the simplest mathematical model 
combining the diffusive (dissipative) mode and the propagative (wave) mode of 
mass transport under local nonequilibrium conditions. In such a case, Eq. (10) 
describes the transport process under non-Fickian diffusion. 

After integration of Eqs. (9) and (10) over an infinitesimal zone that includes 
the interface, the following boundary conditions for the diffusion transport 
hold 

-K L W n T L + K s V n T s = QV n , (11) 

-DV n C = (C - C s )V n + r D -^((C - C s )V n ), (12) 

where Q is the latent heat of solidification, V n Tj and V„C are the normal gra- 
dients of temperature and solute concentration to the interface, respectively, 
V n is the normal velocity of the interface, Cs is the solute concentration at 
the interface in the solid phase given by expression 

C s = kC, (13) 

and k is the coefficient of solute partitioning at the interface. 



2.2 Governing equations 



Under the assumptions drawn in Section 2.1 we consider the governing equa- 
tions for analysis of the morphological stability of the planar interface against 
small perturbations of its form. Our analysis is based on the analysis given 
by Trivedi and Kurz [13] which is advancing the treatment of Mullins and 
Sekerka [10] to the case of rapid solidification. 

The analysis of interface stability is given further in the reference frame mov- 
ing with the constant velocity V. We shall consider the case when the planar 
interface given by equation z(x) = moves along the z-axis of the Carte- 
sian coordinate system (x,z). Then the 2D steady-state fields of the solute 
concentration and temperature are obtained from Eqs. (9) and (10) as 

- 0, (14) 
w. (15) 

a L oz 
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d 2 T s d 2 T s V dT s () 
dx 2 dz 2 as dz 

Following the standard procedure of the analysis of the morphological stability, 
let's place harmonic perturbation on the planar interface. The perturbation is 
described by 

z = (f)(x,t) — 5(t)sm(ujx), (17) 

where 5 is a small amplitude of perturbation (\5\ C l),w = 2n/\ is the cyclic 
frequency with the wavelength, A. The response functions on the perturbed 
interface <p(x,t), i.e. the temperature, T^, and solute concentration, C^, are 
defined by the following relation 

T = T m + mC + FK, (18) 

where m is the slope of the liquidus line in the kinetic phase diagram (i.e., 
phase diagram for nonequilibrium solidification of a binary system), T is the 
Gibbs-Thomson coefficient (i.e. the capillary parameter defined by the sur- 
face energy of the interface), and K is the mean curvature of the perturbed 
interface. 



For the sake of simplicity of the following analysis, the kinetic term V/fj, (in 
which n is the kinetic coefficient of atomic attachment to the interface) is 
omitted in Eq. (18). This simplification has no influence on the main results 
of the present analysis due to the fact that a constant kinetic coefficient does 
not affect the marginal condition of the front stability [11,30]. It has been 
shown that the kinetic effects influence the interface stability in the case where 
temperature dependence of the kinetic coefficient [30], or when an optimum 
stability conjecture for selection of the interface operating state, is used [31]. 

Considering only the terms of the first order of smallness by the amplitude S 
of perturbation, one can obtain for the curvature of the perturbed interface 
the following expression 

Thus, the system (14)-(19) is the extension of the model for the case of sig- 
nificant deviations from local equilibrium in the solute diffusion field during 
rapid solidification. In the case of the local equilibrium solute diffusion, i.e. as 
Vd — > oo, Eq. (14) describes the approximation of Fick's diffusion which has 
been used in the existing models (see, e.g. Ref. [4]). The system (14)-(19) has 
been used in the description of local nonequilibrium solidification with planar 
and non-planar solid-liquid interfaces [19,20,32,33]. With the definition of the 
functions for solute partitioning k(V) and the liquidus slope m(V), the system 
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of Eqs. (14)-(19) can be also applied to the problem of morphological stability 
of the interfaces in rapid solidification. 



3 Morphological stability 



Within the linear analysis of stability, a solution of equations (14)-(16) on the 
perturbed interface (17) is described by 

C,i, = C + b5(t)sm(ux), (20) 



T 4> = T + a5(t)sm(ujx), (21) 
where T and Cq are the values for the temperature and solute concentration 
for the unperturbed planar interface, respectively, and the parameters b and a 
define the corresponding corrections to the small perturbations on it. Taking 
into account the small magnitude of perturbations on the planar interface, the 
perturbed steady-state solution for solute concentration can also be presented 
in a form proportional to S(t) sin(wx). This leads to the following expression 

C(x, z) = C(z) + F(z)5(t) sin(cux), (22) 

where C(z) is the solute distribution for the planar interface. The function 
F(z) is defined from Eq. (20). The far-field condition takes the value for the 
planar interface, i.e. F(z) — > with z — > oo. 

Substitution of Eq. (22) into solute diffusion equation (14) leads to the follow- 
ing approximations: 

- in the zero order of magnitude by the amplitude of perturbation 5 

- in the first order of magnitude by the amplitude of perturbation 5 



A general solution of Eq. (23) has the following form 

C(z) = d + C 2 exp 



D(i-vVyg). 



(25) 



Solution (25) must be limited at z — > oo and it satisfies the following condi- 
tions: C = C and dC/dz = Gc, where Gq is the concentration gra- 

2 = 2 = 

dient at the unperturbed interface. Solution (25) satisfies these conditions in 
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the following form 

G C D{1 - V 2 /V 2 ) 



C(z) = 



C + 



V 



1 — exp 



Vz 



D{i-vyvE), 



, v<v D , 
v>v D . 

(26) 



A general solution of Eq. (24) has the following form 

F{z) ~ F ° exp y 2D(i-vyv$) 

' V-(v 2 + 4D 2 (1-V 2 /V£)uj 2 ) 1/2 



+ F 1 exp 



V 



2D(1 - V 2 jVl) 



(27) 



Solution (27) must be limited at z — > oo and it takes only the real values for 
any u. In this case, Eq. (27) leads to the following particular solution 



F(z) = 



F exp 



u c z 



1/2 



, v<v D , 



(i - vyv 2 ) 

0, V > V D . 



(28) 



Here F is a constant of integration, and the frequency ojc is related to the 
frequency uj of the perturbation as 



V 



UJ C 



2D(i - vyv^y/ 2 



V 



2D(l-V 2 /Vg)V 2 t 



1/2 



(29) 



Note that once the interface velocity is equal to or greater than the diffusion 
speed, V ^ V D , solution (27) takes the zero values for both constants of inte- 
gration. F = is due to limiting of the solution at the infinite point z — > oo, 
and F\ = is for getting the real solution of F(z). Hence, from Eqs. (26) and 
(28) one can obtain that, with V ^ Vd, the coefficient b = in Eq. (20). 

For obtaining F it is necessary that Eq. (22), after substitution of Eqs. (26) 
and (28), must satisfy to Eq. (20) on the perturbed interface (17) in the first 
order of magnitude. It leads to the following expression 



F = b - G 



c- 



(30) 



Hence, substituting Eqs. (26) and (28)-(30) into Eq. (22), one gets an expres- 
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sion for the perturbed field of solute concentration. This yields 



f G C D(1 - V 2 /V*) 



C-C 



V 



1 — exp 



Vz 



+ {b — Gc)S(t) sin(ux) exp 



D(i-vyv*)j 

■(l_V2/Vg)l/2 



v < V D , 



V > V D . 



(31) 

Within the local equilibrium limit Vd — > oo (i.e. when the interface velocity is 
much smaller than the diffusion speed, V <C Vd), solution (31) transforms into 
solution obtained in Ref. [13] for the case of the local equilibrium solute diffu- 
sion transport. Furthermore, as Eq. (31) shows, with the velocities equal to or 
greater than the solute diffusion speed the concentration does not depend on 
the interfacial perturbations and is equal to those one at the planar interface. 
It is known from solution given in Ref. [19] that the solute concentration at 
the unperturbed planar interface is described by 



( l-k 



C{z) - = 



k 



Coo exp 



Vz 



D(l - vyvE) 



lo, 



v <V D , 



V^V D , 



(32) 



where Coo is the solute concentration in the liquid far from the interface (i.e., 
nominal concentration of the alloy). Then, from Eqs. (31) and (32) it is clear 
that 

C(x, z) = C = with V ^ V D . (33) 

Therefore, even in the presence of perturbations, a transition to the complete 
partitionless solidification, C(x, z) = Coo, proceeds in the alloy with the finite 
velocities V ^ V D . 

A solution of Eqs. (15) and (16) for the heat transfer with the condition (21) 
is obtained in Ref. [13] and has the following form 



Tr - T n 



G L a L 
V 



1 — exp 

a L J 



(34) 



+ (a — GL,)6(t) sin(aj:r) exp(— Ulz), 



rp rp G s a s [ / Vz\] 

Ts - To = ^[ 1 - eXI> {-^)\ (35) 
+ (a — Gs)S(t) sin(ux) exp(usz), 

where Gl and Gs are the gradients of temperature in the liquid and solid on 
the unperturbed planar interface, respectively. The frequencies ul and ojc are 
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described by 



V 

= ~ h 



UJ S 



2a L 

V 
~2a~ s 



+uj 2 

2a L 



+ 



V_ 
2a~, 



1/2 



1/2 



(36) 
(37) 



Using the transport balances (11) and (12), the boundary conditions on the 
perturbed interface <fi{x,t), Eq. (17), are obtained for the steady-state regime 
of solidification as follows 



-Kj 



dTj, 



dz 



dz 



QV, 



-D(l - V 2 /V 2 D ) 



dC_ 

dz 



(l-k)VC 



(38) 



(39) 



Here V is the velocity of the perturbed interface. Substitution of the expres- 
sion V = V + (dS/dt) sin(ux) for the velocity of the perturbed interface into 
Eqs. (38) and (39) gives the condition of stability regarding the sign of the 
function 5~ l d5/dt. The concrete form of the function 5~ l d8/dt is obtained as 
follows. 



In the zero order of magnitude by the amplitude 5{t) of perturbation, substi- 
tution of Eqs. (19)— (21) into Eq. (18) gives the following relation 

T = T m + mC , (40) 

for the temperature To and solute concentration Co on the unperturbed pla- 
nar interface. Eq. (40) is also consistent with the liquidus line in the kinetic 
diagram of phase state. In the first order of magnitude by the amplitude S(t), 
one gets a relation between coefficients a and b in Eqs. (20) and (21). This 
yields 

a = mb- Yu 2 . (41) 
Equation (41) gives a linear relation for the corrections to the temperature 
and solute concentration on the perturbed interface and uses Gibbs-Thomson 
effect for the curved interface. For the heat balance at the interface, let us now 
substitute V = V + (dS/dt) sm(ux) for the interface velocity into the condition 
(38). Then, in the zero order of magnitude by the amplitude 5(t), one gets a 
relation for the temperature gradients and the velocity of the unperturbed 
planar interface. This is 

-K L G L + K S G S = QV. (42) 
In the first order of magnitude one gets the following expression 

Q Ydl = a{KLUJL + KsUs) + K L G L{V/a L - u L ) - K s G s (V/a s + u s ), (43) 
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which gives the change of the amplitude S(t) of perturbation in time according 
to the heat balance. Substituting expression V = V + (dS/dt) sin(cjx) into the 
mass balance (39), one can get in the zero order of magnitude a relation for 
the gradient of concentration, interface velocity and solute concentration in 
the following form 



-D(l - V 2 /V 2 )G C = (1 - k)VC . 



(44) 



According to Eq. (39), the change of the amplitude S(t) of perturbation in time 
is defined by the first order of smallness by the amplitude of perturbation. This 
yields 



(Co - 2DG C V/V 2 ) 



IdS 
5di 



+ G C 



Du c {l - V 2 /V 2 )^ 2 - (1 - k)V 
V - Du c (l - V 2 /V%) 1/2 



(45) 



Note that due to introduction of the finite speed V D of solute diffusion into 
the model, an additional term — 2DGqV /V^ has appeared in the left-hand 
side of Eq. (45) in comparison to the analysis in [13]. The prefactor (Co — 
2DG C V/V D V ) has a positively defined value: one gets (C - 2DG C V/V%) > 
due to negative concentrational gradient Gc ^ with k < 1 at any ve- 
locity V < V D . Then, to obtain limiting cases, one can take the expres- 
sion for Gq from the balance of solute concentration (44) in the form Gc = 
-(l-k)VC /(D(l-V 2 /V%)), and one gets the relation (C - 2DG C V/Vg) = 
C (l + 2(1 - k){V 2 /V%)/{l - V 2 jVl)). From this it follows the first positive 
asymptotic: (C — 2DG C V/V^) = C > with V D — > oo. For obtaining the 
second asymptotic with V — > Vd, let us consider the interface velocity in the 
vicinity of the diffusion speed, i.e. the velocity V = Vd — s, ^ e <C 1. Then, 
for V — > Vd, one can obtain the prefactor in the form 



C - 



2DG C V 



V 2 



Co 1 + 



2(1- k(V))V 2 /V 2 * 



V=V D -s 



vyv 2 



V=V D -6 



= C 1 + (1 



k(v )) Vd +v dk{v) 



2(l-k(V D )) 



v=v D 



From this expression follows two consequences. First, in the case of partition 
solidification with V = V D , one gets k(V D ) ^ 1. In this case, the prefactor 
tends to infinity with e — > 0. Second, with the complete solute trapping, 
&(Vd) — 1 an d dk(V)/dV ^ one gets for V — > Vd that the prefactor has a 
positive sign and limited magnitude for the second asymptotic. This yields 



0<C 1 + 



2(l-k(V))V 2 /V 2 ' 

i - vyv 2 , 



Co l + V 



D~ 



dk(V) 



dV 



< oo. 



V=V D ' 



Consequently, (Co — 2DGcV/V%) is positively defined also for the local equi- 
librium case Vd — > oo and local nonequilibrium case V — > Vd if the complete 
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solute trapping occurs, fc(Vb) = 1, and influences only the speed of decreas- 
ing/increasing of the amplitude of perturbation, but not a selection of the 
stable mode itself. 

From equations (43) and (45), taking into account Eq. (41), the expression 
for the function 5~ 1 dS/dt can be obtained. The sign of 5~ 1 d8/dt defines the 
condition of decreasing, 5~ 1 d5/dt < 0, or increasing, 5~ l d5/dt > 0, of the 
interfacial perturbation in time. With 5~ 1 d5/dt = one has the marginal 
(neutral) stability of the interface [10,13]. 



4 Marginal stability 



4-1 The criterion of marginal stability 



We now consider marginal stability for the neutral stability of a small pertur- 
bation on the planar interface, 5~ l d5/dt = 0. From equations (43) and (45), 
one can obtain expressions for the values of a and b. These are 



K L G L 



lu l - V/a L 
K L u L + K s u s 



+ K S G S 



uj s ~ V/as 
K L u L + K s u s ' 



(46) 



G 



c- 



u c - V/[D(1 - V 2 jVl 



2W21 



v<v D , 



'u c -(i-k)v/[D(i-vyv*y/ 2 Y (47) 
o, v > v D . 

The system of equations (41), (46), and (47) allows one to obtain a relation 
for the constant front velocity V and the frequency a; of a perturbation in the 
steady-state regime by excluding of a and b. This relation can be considered 
as a final form for the condition of marginal stability. 



Let's introduce the following functions of stability 

V/a L 



tc 



K L u L + K s u) S ' 

us + V/ as 
K l uj l + K s ujs ' 

cue ~ V/[D(1 - vvyg) 1 / 2 ] 



V <v, 



Di 



(48) 



(49) 



(50) 



u c -(l-k)V/[D(l-V*/Vgy/ 2 ]' 
0, V ^ V D . 

The functions £x and £s in Eqs. (48) and (49) coincide with those derived by 
Trivedi and Kurz [13]. However, as Eq. (50) shows, the function £ c of concen- 
trational stability differs from the corresponding function derived previously 
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in Ref. [13]. As it follows from Eq. (50), within the local equilibrium limit, 
Vd — > oo, one gets the special case £c = (u c — V/D)/(u>c — (1 — k)V/D) 
which coincides with the result given in Ref. [13]. When V ~ Vd, the func- 
tion £ c given by Eq. (50) takes corrections for the relation of the interface 
velocity, V, and diffusion speed, Vd- With V ^ Vd, the exact equality £c = 
takes place. This equality is the consequence of solution of the problem of local 
nonequilibrium solute diffusion which takes into account the finite speed Vd in 
the bulk liquid. Thus, after substituting Eqs. (46) and (47) into Eq. (41) and 
taking into account Eqs. (48)-(50), one can obtain the criterion of marginal 
stability. This yields 



In the local equilibrium limit, Vd — > oo, criterion (51) transfers into the crite- 
rion of marginal stability obtained in Ref. [13] on the basis of a local equilib- 
rium approach to solute diffusion transport. The introduction of the finite dif- 
fusion speed, Vd, into the model leads to the qualitatively new result, which is 
related to the transition to completely partitionless solidification. As Eq. (51) 
shows, with the finite interface velocity V ^ Vd, the solute diffusion ahead of 
the rapid interface is absent [see solution (31)], and the morphological stability 
is defined by the relation of the stabilizing force Tuj 2 , due to surface energy, 
and the contribution K L G L ^ L + KsGsCs of temperature gradients, Gl and 



Using criterion (51), one can analyze qualitatively two different situations for 
solidification when (i) the latent heat is removed from the interface inside the 
undercooled liquid phase, and (ii) the latent heat is removed from the interface 
through the solid crystal phase. In case (i), one gets K L G L £ L + K s Gs^s < 0, 
and the temperature gradient is destabilizing the interface. Therefore, if the 
absolute morphological stability is not reached by the steady balance Too 2 = 
K l Gl^l + KsGsis, the interface is unstable, and the resulting interface may 
exhibit a cellular-dendritic pattern. In case (ii), with K L G L ^ L + K s G s is > 
0, the total heat flux is directed from the front to the solid phase and the 
temperature gradient, in addition to the surface energy, stabilizes the form of 
the solid-liquid interface. In this case, with V < Vd, the morphological stability 
is related to a destabilizing action of the force mGcCc, directly connected with 
the concent rational gradient G c , and the stabilizing force Tuj 2 + K L G L £ L + 
KsGsis, due to the sum of the total positive temperature gradient and the 
surface energy. With the interface velocity V ^ Vd, solidification leads to the 
chemically partitionless pattern. A destabilizing action on the front is absent 
and the interface itself remains linearly stable against any small interfacial 
perturbation. 




(51) 



G s - 
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4-2 Characteristic size for crystal micro structure 



According to the marginal stability hypothesis suggested in Ref. [34] and de- 
veloped in Refs. [4,14], a characteristic size R selected by crystal microstruc- 
ture in solidification (e.g., the dendrite tip radius) is related to the critical 
wavelength, A, of interface perturbation as 

*-A-£. (52) 

Assuming equality for thermophysical parameters of the liquid and solid, one 
can obtain characteristic size, R, from Eqs. (51) and (52). This yields 
- with V <V D 

R ( r/ g y /2 



with V ^ V D 



R -(-\{G& + G&s)) ■ (54) 



In Eqs. (53) and (54) the following designations are accepted 



1 - " 7772 > (56) 



1 + 



2 



1 + 7 TTv^' (57) 

1 + 



aP 



2 



9k 

Ic = 1 + te, (58) 

with P T = VR/2a and P c = VR/2D as the thermal and solutal Peclet num- 
bers, respectively. As we noted above, Eq. (54) is true only for solidification in 
an undercooled liquid, i.e. when the temperature gradient is negative. For the 
case of solidification in the positive temperature gradient the absolute mor- 
phological stability takes place at the interface velocity, V, smaller than the 
diffusion speed, Vd, in the liquid. 
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5 Absolute stability 



5. 1 Nonisothermal solidification 



Let's consider equation (51) in the limit of large perturbations wavelengths 
A>1, which is true for u <C 1. With this condition, from Eqs. (48) and (50) 
one can yield expansions for the functions £l and £c* in the following form 



2 2 



(u 2 D 2 (l-V 2 /VZ) 



Sl KlV2 , & 



V < V D , 



kV 2 ' (59) 

0, V > V D . 



Substituting these expressions into the criterion of marginal stability (51) at 
Gs = 0, we get the criterion of absolute morphological stability for the planar 
front. This criterion can be written in the form of the following nonlinear 
equation for the velocity, Va, of absolute stability: 

V A = V A r (V) + VZ(V), (60) 

where 

VliV) = -2£ (61) 

is the velocity of absolute thermal stability, and 

v c (v) = £. - Yim^ K Vd (62) 

is the velocity of absolute chemical stability. 

The velocity, Vj, of absolute thermal stability, Eq. (61), shows the relation- 
ship between the contribution of the temperature gradient, Gt, and capillary 
parameter, T. With the negative temperature gradient, Gl < 0, the range of 
morphological stability shrinks. Conversely, the positive temperature gradient, 
Gl > 0, extends the range of the velocities at which the planar front is linearly 
stable. The velocity, VJp, of absolute chemical stability, Eq. (62), defines the 
contribution of the concentrational gradient, Gc, and capillary parameter F. 
This velocity is always less than the diffusion speed, < Vd, because the 
velocity V% is defined by the steady balance between surface tension, given 
by the capillary parameter T, and the gradient Gc of solute concentration, 
existing up to the completion of solute diffusion. Consequently, Eqs. (60)-(62) 
exhibit a competition of destabilizing and stabilizing forces. With the velocity, 
V < Va, the planar interface is perturbed with a possible originating of the 
cellular-dendritic patterns. As the solidification velocity increases, V > Va, 
the planar interface becomes morphologically stable against any small pertur- 
bation of its form. 
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To clarify contributions from both thermal and solute diffusion on the absolute 
stability of the interface, we define the gradients in Eqs. (61)-(62) in explicit 
form. For the thermal and concentrational gradients at the unperturbed planar 
interface we use a solution of the local-nonequilibrium problem [19]. From the 
solution, one gets 



Gi 



TqV 
a L 



G c 



(1 - k)VCoo 

kD(i-v 2 /v%y 



v<v D , 
v^v D , 



(63) 



where Tq is a unit of undercooling equal to Q/xl, and is the solute con- 
centration in the liquid far from the interface. From the second expression 
in Eq. (63), it follows that when the solute diffusion is absent ahead of the 
interface with V > Vp, the gradient of the solute concentration is zero ex- 
actly. Substitution of Eq. (63) into Eqs. (60)-(62) gives the expression for the 
absolute stability of the interface. This yields 



V A 



vl + vz, V<V D , 



rC 



V^V D , 



( ^AT T + ^-AT c , V<V D , 



a L 



V^V d , 



(64) 



where AT? = Tq is the thermal undercooling, which is necessary for solidifica- 
tion with the planar interface on the thermal scale, and ATc — (k — tymC^/ k 
is the constitutional undercooling, which is necessary for solidification with 
the planar interface on the scale of solute diffusion. Additionally, AT C is 
the nonequilibrium temperature interval of solidification between liquidus and 
solidus lines in the kinetic diagram of phase state. 

The criterion Vj = cllATt/T in Eq. (64) is the same as that which has been 
obtained by Trivedi and Kurz [13] using the advanced model for large growth 
velocities. The criterion = DAT C /T = D(k - l)mC 00 /(rA; 2 ) in Eq. (64) is 
similar to that which has been obtained by Mullins and Sekerka [10] for the 
case of small growth velocities, and re-derived by Trivedi and Kurz [13] for 
the case of rapid solidification. In addition to this treatment, by introducing 
the finite speed, Vjj, into the model we reach a qualitative new result. At the 
finite velocity, V > V D , due to the absence of the solute diffusion [G c = 0, 
Eq. (63)], the interval between nonequilibrium liquidus and solidus lines is 
equal to zero, ATc — 0. These lines converge in the kinetic phase diagram 
with V >Vd [19], and the absolute stability of the planar interface is defined 
only by the undercooling AT?, and relation between the thermal diffusivity, 
a,L, and capillary parameter, T, Eq. (64). 
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5.2 Isothermal solidification 



In the analysis of the criterion of marginal stability (51), a special interest is 
given to the case in which the form of the interface is defined by the compe- 
tition between stabilizing force, Too 2 , due to surface energy, and destabilizing 
force, mGc^Ci due to gradient, Gq, of solute concentration. Assuming the zero 
temperature gradient, Gl = 0, in Eq. (60), one can obtain an explicit expres- 
sion for the condition of absolute chemical stability of the interface, Va = V% . 
Using the expression for Gc from (63), one gets 

mD(k - 1)Cqo 

A = fjfe2 < (65) 

The form of this expression coincides with the expression given for the case 
of local equilibrium solute diffusion transport at Vp — > oo and V <^Vd [13]. 
However, a final form of the function Va(Coo) is defined by the functions of 
solute partitioning, k(V), and the slope, m(V), of liquidus line in the kinetic 
phase diagram. The behavior of theses functions is rather different for the 
cases of local equilibrium and local nonequilibrium solute diffusion [19,35]. 



6 Discussion 

In the first part of the discussion, we synthesize our system of equations to 
give the self-consistent model, which is adopting the deviation from local equi- 
librium in the solute diffusion field for all functions. We discuss the solute par- 
titioning functions and the expression for the slope of kinetic liquidus which 
take into account the deviation from local equilibrium, both at the interface 
and in the bulk liquid around the interface. Then, in the second part of the dis- 
cussion, we present a quantitative evaluation of the discrepancies between the 
present model and the model in which the local nonequilibrium is taken only 
at the interface. These are compared with experimental data on the absolute 
stability of the planar interface. 



6.1 Solute partitioning and kinetic liquidus 

The boundary condition for solute diffusive transport can be given on the 
basis of the continuous growth model (CGM) [36]. The CGM gives the solute 
partitioning function at the solid-liquid interface, which, in the dilute solution 
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approximation, is described by [36,37] 

uv\ ke ± V I Vdi 

k{V) = T^v/v^' (66) 

where Vdi is the speed of diffusion at the interface, and k e is the value of the 
equilibrium partition coefficient with V — > 0. One of the deficiencies of the 
function (66) is the difficulty to describe of the complete solute trapping at 
the finite interface velocity, i.e. it predicts k — > 1 only with V — > oo. How- 
ever, as it has been shown in numerous experiments (see, e.g., Refs. [38,39]), a 
transition to partitionless solidification occurs at a finite solidification veloc- 
ity. Furthermore, the molecular dynamic simulation has shown [40] that the 
transition to the complete solute trapping is observed at finite crystal growth 
velocity. Therefore, in addition to Eq. (66), a generalized function for solute 
partitioning, in the case of local nonequilibrium solute diffusion within the 
approximation of a dilute alloy, has been introduced [41]. This yields 



fkeil-VyV^+V/V! 



k(V) = < 



DI 



v <V D , 



1-Vyv* + V/V DI ' ' ^ (67) 

i, v > v D , 



where Vdi is the interfacial diffusion speed with Vdi < Vd [20,41]. In the lo- 
cal equilibrium limit, i.e. when the bulk diffusive speed is infinite, Vd — > oo, 
expression (67) reduces to the function k(V), that takes into account the de- 
viation from local equilibrium at the interface only, Eq. (66). In addition to 
the previous model [36,37], the function k(V) described by Eq. (67) includes 
the deviation from local equilibrium not only at the interface (introducing in- 
terfacial diffusion speed Vdi) , but also in the bulk liquid (introducing diffusive 
speed Vd in bulk). As Eq. (67) shows, the complete solute trapping, k(V) = 1, 
proceeds at V — Vd- 

A thermodynamic approach applied to the solidification of a binary system 
[42] provided two models for the solute trapping with and without solute drag 
[23,43]. These models give a shift from local equilibrium at the interface which 
can be expressed in unified form for the slope, m(V), of kinetic liquidus by 
the following equation 

m(V) = |l -*+[* + (1 - k)5 ] In (J^j j. (68) 

Here 5 = is for the model of solute trapping without solute drag and 5 — 1 
is for the model of solute trapping with solute drag. Introducing Eq. (66) into 
Eq. (68), one obtains the constant liquidus slope, m, (independent of V) only 
with the infinite interface velocity, V — > oo. 

Using the results of the local nonequilibrium thermodynamic analysis [35], one 
arrives to the slope of the liquidus line in the following form 
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l — k + lni — \+(l — k) 




2 v_ 

Vd 



v<v D , 



m(V) = < 



(69) 



m e lnk e 



V > V D . 



< k e 1 



With V < Voi the function m(V) includes the function described by Eq. (68) 
for the solute trapping with solute drag (So = 1) and the additional term 
(1 — k) 2 V/Vo- This term arises from the analysis of the Gibbs free energy, 
taking into account local nonequilibrium solute diffusion around the interface. 
It is necessary to note that the function m(V) described by Eq. (69) plays 
a crucial role for self-consistency of the theory of local nonequilibrium solid- 
ification. This form of the function has been used in a self-consistent model 
for rapid dendritic growth and gave quantitative agreement with experimental 
data on kinetics of alloy solidification [20,32]. In particular, the self-consistent 
dendritic growth model, including Eq. (69), predicts the breakpoint at V — Vd 
with good agreement of data on a number of investigated alloys. Furthermore, 
Eq. (69) gives us the ability to describe a transition from the growth kinetics, 
with solute drag effect at small and moderate solidification velocities (aris- 
ing with the developed solute profile ahead of the interface), to the growth 
kinetics without solute drag at high solidification velocities (with the degen- 
eration of the solute profile ahead the interface) [35,44]. Thus, Eqs. (67) and 
(69) close the system of equations (14)-(19) for the self-consistent analysis of 
morphological stability. 

6.2 Comparison with experimental data 

To discuss the results obtained for the interfacial stability, we now compare 
the model predictions for the absolute stability condition (65) in both cases 
of solidification, namely, with local equilibrium solute diffusion and with local 
nonequilibrium solute diffusion transports. Substituting functions (66)-(69) 
into Eq. (65), we analyze the absolute stability of the planar interface for dif- 
ferent velocities. We stress two important points regarding the choice of the 
expression for the slope of kinetic liquidus given by Eq. (68). First, the result 
on rapid dendritic growth [20,32] gives evidence to the confluence of all model 
predictions at small undercoolings and low growth velocities. Disagreement of 
the kinetic curves begins from the undercooling approximately corresponding 
to the undercooling for the absolute chemical stability at moderate growth 
velocities. Therefore, our present discussion for interfacial stability is limited 
by the moderate interface velocities, when the model's predictions (with or 
without taking local nonequilibrium in bulk liquid) begin to disagree. Sec- 
ond, in this region of velocities, the local nonequilibrium approach to rapid 
solidification gives a similar result with the model, which takes into account 
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the deviation from local equilibrium only at the interface, with the solute- 
drag effect (see the analysis presented in Refs. [35,44]). This fact is due to 
the existence of the developed solute profile ahead of the interface at small 
and moderate velocities, when the solute drag may appear at the interface. 
Consequently, in order to evaluate the disagreement between the model and 
experimental data, we choose the expression for the slope of the kinetic liq- 
uidus, which adopts the solute-drag effect, i.e. it is chosen in the following 
calculations where 5 — 1 for Eq. (68). 



Using parameters of an Al-Fe alloy from Table 1, one can calculate the curve 
for critical concentration, C^V), which gives a threshold for interface insta- 
bility. As it can be seen from Fig. 3, two regions of the interfacial existence 
may occur: the planar interface is absolutely stable below the curves and the 
interface breaks down in the regions above the curves given by the functions 
CooiV). In comparison with the model with the local equilibrium diffusion 
and deviation from local equilibrium at the interface only [Eqs. (65, (66), and 

(68) ], the present model for interface stability with the deviation from local 
equilibrium, both at the interface and in the bulk liquid [Eqs. (65), (67, and 

(69) ], defines a curve CooiV) which is limited by the diffusion speed Vd for 
morphological stability of the interface. This limit exists due to a steady bal- 
ance between the stabilizing capillary force and the destabilizing force defined 
by the concentrational gradient, which still acts on the interface up to the 
finishing of diffusion, i.e. until the point V = Vd- 



For a quantitative comparison of the model predictions we have chosen experi- 
mental results on interface stability during rapid solidification of a Si-Sn alloy 
as presented by Hoglund and Aziz in Ref . [46] . These authors measured a crit- 
ical concentration of Sn for interface breakdown in a steady-state solidification 
after pulsed laser melting. Using the parameters of the Si-Sn alloy from Table 
1, the model predictions for the function C^iV) are compared quantitatively 
with experimental results from Ref. [46]. As is shown in Fig. 4, the present 
model for interface stability with the local nonequilibrium diffusion [Eqs. (65), 
(67], and (69)] gives a satisfactory comparison with the experiment. At the 
concentration Sn = 0.02 atomic fraction (see the extreme right experimental 
point in Fig. 4), the discrepancy between the model with local equilibrium 
solute diffusion (curve 1 in Fig. 4) and experiment gives the value of 38.90 % 
(see Table 2). At the same alloy's concentration, the present model (curve 2 in 
Fig. 4) gives the discrepancy with experiment of 16.93 % (see Table 2). Conse- 
quently, even better comparison with the available experimental data can be 
obtained with using the present model of local nonequilibrium solidification. 
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7 Conclusions 



Morphological stability of the planar interface in rapid solidification of non- 
isothermal binary system has been considered. We have taken into account the 
fact that the high rate of solidification process leads to the absence of a local 
thermodynamic equilibrium in the solute diffusion field and at the solid-liquid 
interface. The presently developed model is self-consistent: the main governing 
equations, Sec. 2, and the interface conditions for solute trapping and kinetic 
liquidus, Sec. 6, are consistent with the formalism of extended thermodynamic 
approach to rapid solidification [35] . Using the model of local nonequilibrium 
rapid solidification, our analysis of morphological stability extends the previ- 
ous analysis of Trivedi and Kurz [13], which has been performed to advance 
the treatment of Mullins and Sekerka [10] to the case of rapid solidification. 
The main outcomes of this analysis are summarized as follows. 

(i) For the velocities equal to or greater than the diffusion speed, V > Vd, 
from solutions (31)-(33) it follows that the field of concentration does not 
depend on a form of the interfacial perturbation and it is equal to the initial 
(nominal) concentration, C(x, z) = C^. This result is in agreement with the 
previous results for planar and parabolic interfaces [19,33]. Solutions (31)- 
(33) has a clear physical meaning: a source of concentrational disturbances, 
i.e. the perturbed interface, cannot disturb a binary liquid ahead of itself if 
the interface velocity is equal to or greater than the maximum speed of these 
disturbances. 

(ii) The obtained criterion of the marginal stability, Eq. (51), defines a wave- 
length of perturbation for the neutral stability. For V < Vd, the neutral sta- 
bility is defined by a balance of the stabilizing force, due to surface energy, 
destabilizing force, due to concentrational gradient Gc, and the contribution 
of temperature gradients Gl and G$- Qualitatively, a new result can be ob- 
tained from the criterion for the front velocity of V > Vd, i.e. for an absence 
of the solute diffusion ahead of the interface. As it follows from the second 
equation in Eq. (51), the morphological stability of the interface is defined 
only by the relation between the thermodynamic stabilizing force, due to the 
surface tension, and the driving force of the morphological instability, due to 
the negative thermal gradient in the undercooled liquid. In the case of direc- 
tional solidification with the positive thermal gradient, the destabilizing action 
on the interface is absent at V > Vd, and the interface remains linearly stable 
against small perturbations of its form. 

(iii) Absolute stability of the planar interface is considered as a steady balance 
between destabilizing force (due to the concentrational gradient), the thermal 
contribution (due to the thermal gradient), and the stabilizing force (due to 
surface tension). The velocity, V A , of the absolute interface stability is obtained 
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as a sum of the velocity, Vjf , for thermal stability and velocity, Vj( , for chemical 
stability defined by Eq. (64). is the same as what was obtained by Mullins 
and Sekerka [10] for the case of small growth velocities. It was re-derived by 
Trivedi and Kurz [13] for the case of rapid solidification. Introduction of the 
finite speed, Vb, into the model gives the qualitatively new result: with the 
absence of the solute diffusion at V > Vd, the absolute stability of a planar 
interface is defined only by the thermal undercooling and relation between the 
thermal diffusivity, a L , and capillary parameter, T, Eq. (64). For an isothermal 
solidification, the present analysis shows the limiting boundary equals to the 
diffusion speed, V D) for the region of morphological instability [see Eq. (65) 
and Fig. 3]. 

(iv) The predictions of the present model for the critical concentration above 
which a planar interface becomes unstable [see Eqs. (65), (67), and (69)] are 
compared with the previous model, which adopts the deviation from local 
equilibrium at the interface only [see Eqs. (65), (66), and (68)], and with the 
experimental data obtained for solidification of the Si-Sn alloy [46]. As it is 
shown in Fig. 4, the present model is able to describe experimental data 
satisfactorily in a whole region of the interface velocities investigated. From 
numerical evaluation of the theoretical predictions summarized in Table 2, 
it follows that better agreement with experiment is obtained with using the 
present model of local nonequilibrium solidification. 
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Table 1 

Physical parameters used in calculations of the limit of absolute stability of the 
planar interface in solidifying binary alloys 



Parameter Al-Fe Si-Sn 



Diffusion coefficient 


D L 


m 2 /s 


1.7 x i(r 9 (*) 


2.5 x 


io- 8 


(**) 


Partition coefficient 


k e 




0.03 (*) 


0.016 


(**) 




Liquidus slope 


m e 


K/at.% 


-7.3 (*) 


-4.6 


(**) 




Gibbs-Thomson coefficient 


r 


K • m 


1 x 10~ 7 (*) 


1.3 x 


io- 7 


(**) 


Interface diffusion speed 


Vdi 


m/s 


7 


17 (**) 




Diffusion speed in bulk liquid 


V d 


m/s 


10 


17.5 







(*) - Data taken from Ref. [45]. 
(**) - Data taken from Ref. [46]. 



Table 2 

Discrepancy between theoretical predictions and experiment for the absolute chem- 
ical stability of the planar interface in the Si-0.02(at.fraction)Sn alloy 



Absolute chemical stability 


Velocity 


Definition 


Value 




(m/s) 


of discrepancy 


of discrepancy (%) 


Local equilibrium solute diffusion, 
Eqs. (65), (66), and (68) 


= 15.5 


y-(l) v (exp) 

100% 


38.90 



Local nonequilibrium solute diffusion, 

y(2) _ y(exp) 

Eqs. (65), (67), and (69) V\ 2 ' = 11 A ^ 100% 16.93 



Experiment, Ref. [46] V£ xp> = 9.47 
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Figure captions to the article "Linear morphological stability analysis of the 
solid-liquid interface in rapid solidification of a binary system" . 

Fig. 1. Morphological diagram for solidification of binary systems which is illus- 
trating the microstructural transitions "planar front" — > "cellular structure" — > 
"dendrites" — > "cellular structure" — > "planar front", with the increasing of the 
solidification velocity, V. Here Vc is the velocity given by the criterion of constitu- 
tional undercooling, and Va is the velocity for absolute morphological stability of 
the interface. 

Fig. 2. Optical micrograph of longitudinal through-thickness section of a melt-spun 
ribbon of Ni-18 at.% B [24]. The crystal microstructure exhibits a transition from 
planar interface with solute segregation-free to cellular-dendritic patterns. The tran- 
sition proceeds due to decreasing of the interface velocity from V > Va up to 
V < Va- Wheel surface at bottom of micrograph. 

Fig. 3. Velocity Va of absolute chemical stability versus solute concentration Co 
for Al-Fe alloy. Dashed curve corresponds to solution of Eqs. (65), (66), and (68) 
with solute-drag effect. Solid curve corresponds to solution of Eqs. (65), (67), and 
(69). Dashed-dotted line, V = Vd, represents the limiting velocity for the absolute 
interface stability. 

Fig. 4. Critical concentration, Coo, above which planar interface is unstable. Exper- 
imental points correspond to solidification of the Si-Sn alloy [46]. Circles are taken 
from measurements performed on bulk single crystal Si(100), and squares are taken 
from measurements using Sn-implanted Si-on-sapphire (SOS) samples. Curves are 
given by the models for interfacial absolute stability: dashed - with local equilib- 
rium diffusion and solute drag effect, Eqs. (65), (66), and (68); solid - with the local 
nonequilibrium diffusion, Eqs. (65), (67), and (69). Dashed-dotted line, V = Vd, 
represents the limiting velocity for the absolute interface stability. 
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Figure 4 
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